INP MSU 98-13/514 
March 23, 1998 



Higgs Boson Discovery Potential of LHC 
in the Channel pp ^ 'yy -\- jet 



S. Abdullin+, M. Dubinin*, V. Ilyin*, 
D. Kovalenko*, V. Savrin*, N. Stepanov'^ 

+ UHA, Mulhouse, France (on leave of absence ITEP, Moscow, Russia) 
Institute of Nuclear Physics, Moscow State University, 119899 Moscow, Russia 
^ Institute for Theoretical and Experimental Physics, 117259 Moscow, Russia 



Abstract 



We discuss the SM Higgs discovery potential of LHC in the reaction pp 
H + jet — >■ 77 + jet when the jet is observed at sufficiently high Et to be reliably 
identified. We conclude that this channel gives promising discovery possibilities for 
the Higgs boson mass range 100-140 GeV, during LHC operation at a low luminosity. 
With 30 fb~^ of accumulated data and for = 120 GeV about 100 signal events 
could be observed with the number of background events larger by a factor of 2 only, 
showing a signal significance S/ ^/B ~ 7. We use the difference of distributions in the 
partonic subprocess energy \/§ for the signal and background for a better separation 
of the signal. 



Introduction 



It is well known that the observation of Higgs boson with mass Mjj < 140 GeV at the LHC 
collider (pp, y/s =14 TeV) in the inclusive channel — > 77 + X is not easy [0, |^. The 77 
continuum rapidly increases for smaller 77 pair invariant masses, and it is necessary to separate 
a rather elusive Higgs boson signal from it. In this situation it is important to understand 
whether we can observe any other discovery channels. In this paper we are considering (in the 
Standard Model) the reaction pp ^ H + jet — *■ 77 + jet when Higgs boson is produced with 
large transverse momentum recoiling against a hard jet. Of course, in this channel the signal 
rate is much smaller in comparison with the inclusive pp —>■ '-ff + X case. At the same time, 
as we shall see below, the situation with the background is undoubtedly much better. We can 
usefully exploit richer kinematical features of the final state 77 + jet, when some specific jet 
distributions in the partonic c.m.s. are different for signal and background processes. 

The idea to look for Higgs boson associated with a high Et jet in the final state was con- 
sidered in 1^. In IQ, 1^ the corresponding subprocesses were calculated within the minimal 
supersymmetric extension of the Standard Model. However, in ^ Higgs decay channels were 
not considered and in only the final state t~^t~ + jet was analyzed. In the SM heavy 
Higgs boson decay into the WW or ZZ pairs was discussed for H + jet production. Promising 
results have been obtained recently in for 77 + 2jets final state with two very forward jets 
(|?7| < 5). In P the final state 77 + (> 2jets) was simulated in the realistic CMS detector 
environment. The final state 77 + jet with only one high Et jet recoiling against the Higgs 
boson has not been analyzed yet and we investigate it in detail Q 

1. Calculation framework 

We calculated cross sections and distributions by means of CompHEP package ||T0|. Methods 



of multichannel adaptive Monte Carlo integration over the phase space implemented in the 



package are described in [p]] . 



The cross sections of signal subprocesses under discussion depend significantly on the choice 
of QCD parameters. There are three sources of this dependence: (1) QCD evolution of the 
parton distribution functions; (2) dependence in the subprocesses; (3) if — > 77 branching 
variation due to QCD corrections to the F^*. In the leading order the corresponding corrections 
can be factorized. Moreover, due to a very small value of F^* for < 140 GeV the fixed value 
of strong coupling asiMn) can be used for the evaluation of ^ 77 branching. However, 
it is well-known that for the reaction — ii — 77 the dependence on the renormalization 
scale /i and on the parton factorization scale Q is strong enough, and the next-to-leading 



order analysis is necessary (see ||T2| and references therein). The NLO corrections decrease 
this (yU, Q) theoretical uncertainty, showing only a ~ 15% sensitivity of the final result. One 
can hope to observe a similar effect also for the case of Higgs boson production at high pt. 
However, self-consistent analysis requires the NLO corrections to hard subprocesses which are 
not known yet, unfortunately^. Thus, today we made only the LO analysis when parton 
distribution functions and running in subprocesses are taken at the leading order. We 
used the parametrization CTEQ41 (^.(M^) = 0.1317, AgJ.^ = 181 MeV) At the same 

time QCD NLO formulas |jT5|, |12|] were used for the evaluation of F^* with the reference value 
Br{H 77) = 1.534 ■ 10"^ at Mh = 100 GeV. Parameter values used in our analysis are 

^Primary analysis of the reaction pp 77 + jet can be found in where we also considered the processes 
pp — > WH, UH ^ 77 + lepton. 

^Recently formulas for the corresponding virtual corrections were calculated in the mt 00 limit 
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Mz = 91.187 GeV, sin^^ = 0.4732, and m, = 0.2 GeV, = 1.42 GeV, = 4.62 GeV, 
mt = 175 GeV. 

Typical acceptances of the LHC detectors ATLAS and CMS should be taken into account 
in the analysis. For photons we are using the cuts similar to already well analysed cuts of the 
inclusive channel [|I], 0] : two photons are required with p] > 40 GeV for each photon (somewhat 
harder than for the inclusive channel), and photon rapidity \r]^\ < 2.5. For a jet two sets of 
basic kinematical cuts will be discussed in further analysis: 

CI: > 40 GeV, \r]jet\ < 2.4; C2: > 30 GeV, \7]jet\ < 4.5. 

The first set corresponds to the central part of calorimeter where a jet with transverse energy 
greater than 40 GeV can be detected with highest efficiency. The second set assumes forward 
parts of the hadron calorimeter involved. The analysis of background processes made for the 
inclusive channel has shown that photons should be isolated. So we apply the cut AR > 0.3 
for each 7 — 7 and 7 — g{q) pair, where AR = a/A0^ + Ai]'^ is the standard variable separating 
particles in the azimuth angle - rapidity plane. 

2. Signal processes 

There are three QCD subprocesses giving a signal from the Higgs boson in the channel 
under discussion: gg —>■ H + g, gq ^ H + q and qq —>■ H + g. Feynman diagrams contributing 
in the leading order are shown in Fig. [l|. The corresponding matrix elements were calculated 
analytically in and we implemented these formulas in our code. We found that the gg 
H + g "yj + g subprocess gives the main contribution to the signal rate. The gq channels with 
light quarks give 20-25% of the gluon-gluon signal contribution, while the qq channels can be 
neglected. Furthermore, at large values of parton factorization scale, ~ 10^ GeV^, typical for 
intermediate Higgs boson production, one can expect noticeable contributions of the strange 
and heavy quark sea. We found that it gives 10% of the dominant gluon-gluon subprocess. 
Then, tree-level subprocess gg —>■ H + b -\-b, where the Higgs boson is radiated from the final 
b quark or ^-propagator, contributes to 77 + jet signal events if we do not register {veto) the 
events with both b quark jets observed. We found that this process gives about 1.5% of the gg 
contribution. Finally we note that the s-channel subprocesses qq ^ H + Q + Q, where g is a 
light quark and Q is 6 or t quark, are negligible. In total QCD subprocesses give 3.3 fb, 5.7 fb 
and 5.5 fb with the CI set of cuts for = 100, 120 and 140 GeV, correspondingly, and 5.5 
fb, 10.6 fb and 9.8 fb with the C2 set. 

The second group of signal subprocesses includes the electroweak reactions of Higgs pro- 
duction through WW or ZZ fusion (Fig. 0a) with the high Et jet in the final state from the 
scattered quark, and the Higgs boson production associated with W or Z (Fig. ||b) decaying 
into quark-antiquark pairs. Of course, if the signature with only one jet is under discussion. 



one should veto the second quark jet. The WW/ZZ fusion processes were calculated in [ p!6[| . 



while the processes of HW/HZ associated production were considered in WTl. One should note 



that the contribution of these subprocesses decreases when we change CI set of cuts to C2, 
what is opposite to the case of QCD signal subprocesses. The reason is veto condition for one 
of the jets. Indeed, typical transverse momentum of the final quark (in the dominant WW/ZZ 
fusion production processes) is about half of the vector boson mass. So, weaker jet cut Et > 30 
GeV in the veto condition removes a larger number of events. We found that the cross section 
of fusion processes is 2-3 times larger than for the processes of associated production. In EW 
channels s and c quark sea gives about 10% of u and d quark (valence plus sea) contribution. 
In total, EW processes give 1.0 fb, 1.6 fb and 1.5 fb with the CI set of cuts for Mh = 100, 
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120 and 140 GeV, correspondingly, and 0.65 fb, 1.1 fb and 1.1 fb with the C2 set. So the 
electroweak signal rate is on the level of 30% of the QCD signal with the CI set of cuts and 
10% with the C2 set. 



3. Irreducible background processes 

Let us now look at the irreducible background. Three subprocesses contribute here, two of 
them at the tree level, gq —>■ 'J + 'y + q (Fig. 0a) and 5'^ — > 7 + 7 + (? (Fig. 0b), and the third one 
at the one loop level, gg —>■ J + J + g- The qq channel with u and d quarks in the initial state 
gives the cross section on the level 15-20% of the main irreducible background gq {q = u, d) 
with the CI set of cuts, and 15-30% with the C2 set for Higgs mass range 100-140 GeV. The 
contributions of strange and heavy quark sea in the initial state are on the level 25-30% of the 
main gq channel. Here about 85% of the contribution is coming from the c quark sea. The s 
quark contribution is suppressed by the factor 16 originating from a smaller fractional charge. 

The matrix element for the one-loop process gg 7 + 7 + ^7, when the photons are radiated 
from the quark loop, is unknown yet. We estimated this cross section by means of PYTHIA 



18| simulation, switching on the gluon bremsstrahlung from the initial state in the subprocess 



gg •y + •y. This is definitely only one of physical contributions to the one-loop process. The 
result of these simulations shows us that the one-loop background is about 2-4% of the main 
contribution coming from the gq channel. So we neglect the one-loop background in further 
analysis, however, this point is one of a serious theoretical uncertainties. 

In total, the irreducible background contribution amounts to 9.3 fb, 13.7 fb and 16.1 fb in 
the 1 GeV bin around the values of M^^ = 100, 120 and 140 GeV, correspondingly, and with the 
CI set of cuts. Of course, with the C2 set of cuts the background is larger, the corresponding 
numbers are 15.3 fb, 25.3 fb and 25.1 fb in the 1 GeV bin. 



4. Reducible background processes 

Various processes could give a background due to radiation of photons from the fragmentating 
quarks or gluons. Photon production in a jet hadronization is also possible, and is defined in 
particular by a vr^-meson production. The energetic vr^'s, decaying mainly to a photon pair, 
can be detected as a single photon in the electromagnetic calorimeter. First kind of such 
reducible background is coming from the subprocesses gq^y + g + q, gg^y + q + q and 
qq' ^ 1 + q{g) + q'{g), in the case when the final gluon or quark produces an energetic isolated 
photon but other products of hadronization escape the detection as a jet. One can say that 
a jet is misidentified as a photon. Second kind of reducible background could come from the 
subprocesses gq j + q, qq 1 + g when the second photon is produced in the quark or gluon 
fragmentation but other products of the hadronization are still detected as a jet with proper 
separation from this photon. Third source of reducible background is connected with the pure 
QCD subprocesses of 2 — 2 and 2 — 3 types, when all particles in the final state are gluons 
or quarks. We performed rough analysis of these QCD reducible backgrounds in and found 
that they turn out to be less than 20% of the irreducible background, and the misidentification 
rate is given mainly by the processes of the first kind. Here we used the 'y{'n'^)/jet rejection 
factor equal to 2500 for a jet misidentified as a photon and 5000 for a well separated 7(vr°) 
production by a jet. These factors were obtained for a jet satisfying the cuts described above 
with the help of PYTHIA simulations. No additional vr*^ rejection algorithms were used. 

5. Dependence on the scale 
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In our calculations we used = Mfj + 2{El^^Y as the parton factorization scale and the 
normalization scale for running as in the hard QCD signal subprocesses, as well as in the 
background processes. At the same time for the WW/ZZ fusion we used = (Mv/2)^ 
(My = Mw, Mz) and = {My + Mnf for the HW/HZ associated production. Certainly, 
due to order of partonic subprocesses one can expect strong dependence on the choice of 

for the QCD signal. We checked that for = Mfj, the corresponding signal cross section 
increases by 15%, while the background cross section increases only by ~ 5%. If one uses 

= (50GeV)^, the QCD signal cross section increases by 80% while the background only by 
13%. Such strong Q"^ dependence shows that the complete NLO analysis is needed. 

6. Reconstruction of event kinematics 

The distributions da/d^/§ in parton collision energy presented in Fig. ^ show that the back- 
ground processes contribute at a smaller in comparison with the QCD signal processes. So, 
one can hope that the corresponding cut can improve the S/B ratio. Photon energy can be mea- 
sured with a high enough presision, so the main uncertainty of the reconstruction is defined 
by the accuracy of a final parton energy and momentum reconstruction (see details in [0, Q). 
The energy of a parton can be reconstructed with the accuracy {5E/EY ~ 100% a/E © 5%. 
It means that for jets with transverse energy more than 30 GeV the parton energy will be 
reconstructed with the accuracy ~ 7 GeV. Taking into account all factors, one can expect that 
\/I variable can be reconstructed with the error ~ 10 GeV. It is clear from the distributions 
presented in Fig. | that such accuracy is good enough to apply a cut on for the suppression 
of the background. However 5E is mostly related to the energy loss due to the undetected 
products of hadron fragmentation. It follows that the experimental distributions can be shifted 
to the smaller values in comparison with Fig. ^. But the scale of smearing in this distribu- 
tion (more dangerous for our analysis) should be definitely smaller than 10 GeV. One should 
note that processes at higher order (contributing to the NLO corrections) have more than 3 
particles in the partonic final state. So, strictly speaking, the energies and momenta of photons 
and the detected jet are not sufficient for the reconstruction of partonic c.m.s. collision energy 
in these events. However, one can hope that the main contribution to the QCD corrections 
should go from virtual and soft gluons and cannot affect significantly the reconstruction of 
variable. 

We found that the cut > 210 GeV decreases the QCD signal cross section only by 25-30% 
with the CI set of basic cuts, while the background is suppressed much stronger, for example 
by a factor of 2 when = 120 GeV. For the C2 set of basic cuts the corresponding numbers 
are 20-35% decrease for the signal rate and suppression factor of 1.7-2.8 for the background. 
It means, for example, that with the C2 set of cuts the S/B ratio is improved by a factor of 
~ 2 for Mh = 100 — 140 GeV. One can get S/B even better by applying much stronger s cut. 
For example \/S > 300 GeV in the case of C2 set suppresses the background by a factor of 8.7 
while the QCD signal only by a factor of 2.6. 



Mh = 120 GeV 


no V§ cut 


> 210 GeV 


VI > 300 GeV 


aQ'^^fb 


10.6 


7.0 


4.0 


af^fb 


1.1 


0.81 


0.58 


ae fb/GeV 


25.3 


9.1 


2.9 



The reducible QCD background should be also suppressed by the a/I cut, even stronger than 



4 



the irreducible one. Indeed, the subprocesses of reducible background have a similar diagram 
topology. So our arguments (see next section) should work in this case also. However, a shift 
to smaller for the reducible background processes will be larger than for irreducible ones, 
because the energy of parton misidentified as a photon will be reconstructed with a much higher 
loss of energy in comparison with the photon radiated from the hard subprocess directly. 

7. Angular distributions in the partonic c.m.s. 

Useful information is provided by the spin structure of in and out-states. For the dominant 
subprocess gg H + g a. set of possible in spin states does not include spin 1, while the 
spin of the out- state is determined by the gluon. It means, in particular, that the S-wave 
does not contribute here. At the same time for the background subprocesses gq ^ 'J + ■y + q 
and qq ^ 1 + 1 + g the same spin configurations are possible for both in and out states. So, 
the jet angular distribution in the partonic c.m.s. should be different for the signal and the 
background. 

In Fig. ^ we represent the angular distributions for the signal and background processes. 
Here d*^^ is the jet scattering angle in the partonic c.m.s. The background curve has a bump 
while the signal curve is smooth enough. In Figs. ^ and ^ the same distributions are shown 
when the cut is applied. One can see that events from the central part are suppressed. 
Qualitative interpretation of this effect is given by the simple observation that the relative 
contribution of partial waves with higher angular momentum increases with the collision energy 
for processes with t-channel exchange of virtual particles. Thus, the contribution of S-wave in 
the background is relatively suppressed, substantiating the effect of the cut on partonic collision 
energy \/I discussed in the previous section. 

It is clear that the angular distributions of photons from the signal and background processes 
in the partonic c.m.s. should be different. Indeed, photons from the Higgs decay (giving uniform 
angular distribution in the Higgs rest frame) are produced mostly in the direction opposite to 
the jet in the partonic c.m.s. At the same time photons radiated from the final quark in gq 
subprocesses, which give dominant contribution to the background, will be observed mainly at 
small angles with the jet. This effect is clearly seen in Fig. ^a, where the distributions in "i?*^ 
(angle between the jet and the photon with smaller pt) are represented in partonic c.m.s. Let 
us now apply the cut. The distributions in this case are represented in Figs. |]b and ||c. 
One can see that events at small angles ■(9*^ are suppressed affecting mainly the background 
distribution. This is qualitatively understood because events with small angle ■(9*^ correspond 
in average to smaller momentum of the Higgs boson or 77 system. Indeed, for high enough 
velocity of the 77 system the corresponding Lorentz boost will turn (almost) all photons in the 
directions opposite to the jet momentum. So, events where one photon is radiated in the jet 
hemisphere should have smaller \/^ than events where both photons are radiated opposite to 
the jet. In particular, this effect gives us one more argument why the a/J cut suppresses the 
background stronger than the signal. 

8. Estimates for the LHC detectors 

In this section we give some estimates of the Signal/Backgroundiatio and the signal significance 
for LHC detectors in the channel — > 77 + jet using basic numbers from ATLAS and CMS 
Technical Design Reports 0. First, the efficiency of photon detection should be taken into 
account. This important instrumental parameter is expected to be ~ 80% for both detectors, 
and it was obtained from the simulation of pp ^ if — 77 events in the detectors. Second, the 
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fiducial area cuts should be considered to exclude the regions of electromagnetic calorimeter 
where the performance decreases crucially. The corresponding efficiency per photon can be 
taken at the level 95% for ATLAS and 92.5% for CMS. Finally, let us consider the M-^-^ bin 
optimization which affects rather strongly the signal significance. At LHC low luminosity 
regime and for = 100 GeV the 80% signal events reconstruction efficiency corresponds to 
the M^^ mass bin of 3.1 GeV (resolution parameter am = 1-1 GeV) for ATLAS Lead-Liquid- 
Argon electromagnetic calorimeter [Q. For CMS PbW04 electromagnetic calorimeter the 
mass resolution parameter am = 0.65 GeV assumes the mass bin AM^^ = 1.9 GeV with the 
reconstruction efficiency 73% in this window. For larger Mfj the M^^ resolution is slightly worse. 
Let us take the values AM^^ = 3.25 GeV and 2.0 GeV if Mh = 120 GeV, and AM^^ = 3.4 
GeV and 2.1 GeV if Mh = 140 GeV, for ATLAS and CMS M^^ bin, correspondingly. 

Then, the QCD next-to-leading corrections should be taken into account. As we have 
mentioned above, these corrections are unknown for the processes discussed here. Some estimate 
may be possible since it was shown (see [jl2[ and the references therein) that the NLO corrections 
to the gg fusion subprocesses in the pp — > 77-l-X inclusive channel increase the Higgs production 
cross section by a factor K ~ 1.6. One can hope that the corresponding iiT-factor in our case will 
be at least of the same value. Of course, it is very probable that the NLO corrections enhance 
somehow the rates of background processes as well. So, in the absence of theoretical results let 
us use in our case the factor K^^^ = 1.6 both for the signal and background subprocesses. 

Finally we get the following numbers in the case of C2 basic kinematical cuts and applying 
additional cut y/I > 300 GeV: 



Low 




100 GeV 




120 GeV 


Mh = 140 GeV 


luminosity 


ATLAS 


CMS 


ATLAS 


CMS 


ATLAS CMS 


al" fb 


1.8 


1.6 


3.4 


2.9 


3.6 3.1 


a%" fb 


6.6 


3.8 


8.9 


5.2 


12.3 7.2 


S/B 


1/3.7 


1/2.5 


1/2.7 


1/1.8 


1/3.4 1/2.3 


S/VB (30 fb-i) 


3.8 


4.3 


6.2 


7.0 


5.6 6.3 



9. Conclusions 

The channel 77 + jet with the jet transverse energy Ef > 30GeV and rapidity \t]\ < 4.5 
(thus involving forward hadron calorimeters) gives very promising discovery possibilities for 
the Higgs boson with a mass of 100-140 GeV during the LHC operation at a low luminosity of 
~ 10'^'^cm~^s~^. For example, with an integrated luminosity of 30 fb~^ about 100 signal events 
could be observed for Mh = 120 GeV with a number of background events only 3 times higher in 
ATLAS and 2 times higher in CMS detector. These numbers demonstrate the main advantage 
of this channel in comparison with the inclusive reaction pp —>■ 77 + X, namely significant 
improvement of the S/B ratio. The estimate of the reducible background using only isolation 
criteria shows that it is less than 20% of the irreducible background. Our results for the signal 
and background rates mean that the discovery level S/ ^/B = 5 for the signal significance will 
be achieved already with an integrated luminosity of 30 fb~^ for Mh = 110 — 140 GeV. 

We found that the detection of jets with rapidities up to \ri\ = 4.5 improves noticeably the 
signal significance. For example, 20-25% improvement can be achieved in comparison with the 
case when hard jets are only centrally produced {\ri\ < 2.4, Ef'^^ > 40 GeV). 
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We demonstrated that sufficiently rich kinematics of the three particle final state allows to 
introduce new observable distributions, suitable for a better separation of the signal. Further 
improvement of the signal significance can be achieved by using the jet and jet-photon angular 
distributions in the reconstructed partonic c.m.s. These distributions can be also used for a 
further suppression of the reducible QCD background. 

One should also note that in the present analysis we used parameters obtained from the 
simulation of the reaction ^ 77 + X in ATLAS and CMS detectors. However, the 77 + jet 
kinematics is more preferable for the event reconstruction than in the inclusive case. Recon- 
struction of the jet in the hadronic calorimeter allows to determine more precisely the position 
of interaction vertex. Photons from the Higgs decay in 77 + jet state arc more energetic than 
for the inclusive channel. So the photon reconstruction efficiency and effective mass resolution 
used in our analysis probably are too pessimistic. Prom the experimental data processing point 
of view, additional event selection criterion (trigger condition) of a jet in the final state allows 
to restrict the number of diphoton events in comparison with the inclusive channel, providing 
opportunities for a more careful analysis at a better S/B ratio. AUtogether these factors could 
give a sizeable improvement of the signal significance. 
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Figure captions 



Figure 1: Feynman diagrams for the QCD signal subprocesses: a) gg ^ H + g , h) gq ^ H + q 
and c) qq ^ H + g. 

Figure 2: Feynman diagrams for the electroweak signal subprocesses with a) WW and ZZ 
fusion mechanisms of the Higgs boson production, b) associated HW and HZ production. 

Figure 3: Feynman diagrams for the subprocesses a) gig — > 7 + 7 + g and b) qq ^ ■y + J + g 
contributing to the irreducible background. 

Figure 4: Distributions in the parton cm. energy v^I for the signal (S) (without EW contri- 
bution) and background (B) processes. Here Mh = 120 GeV and the basic set of kinematical 
cuts C2 is imposed. The 77 invariant mass for the background is integrated over the 1 GeV 
bin. 

Figure 5: Distributions in the jet angle in partonic c.m.s. for the signal (S) (without EW 
contribution) and background (B) processes in the case Mh — 120 GeV. Upper plot is obtained 
with the basic set of kinematical cuts C2. Next plots show changes in this distribution when 
the \/§ cut is applied. The 77 invariant mass for the background is integrated over the 1 GeV 
bin. 

Figure 6: Distributions in the angle between jet and the photon with smaller transverse momen- 
tum in partonic c.m.s. for signal (S) (without EW contribution) and background (B) processes 
in the case Mh = 120 GcV. Upper plot is obtained with the basic set of kinematical cuts C2. 
Next plots show changes in this distribution when the \/§ cut is applied. The 77 invariant 
mass for the background is integrated over the 1 GeV bin. 
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